Vibrational wave packet induced oscillations in two-dimensional electronic spectra. 

II. Theory 
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We present a theory of vibrational modulation of two-dimensional coherent Fourier transformed 
electronic spectra. Based on an expansion of the system's energy gap correlation function in terms 
of Huang-Rhys factors, we explain the time-dependent oscillatory behavior of the absorptive and 
dispersive parts of two-dimensional spectra of a two-level electronic system, weakly coupled to 
intramolecular vibrational modes. The theory predicts oscillations in the relative amplitudes of 
the rephasing and non-rephasing parts of the two-dimensional spectra, and enables to analyze time 
dependent two-dimensional spectra in terms of simple elementary components whose line-shapes are 
dictated by the interaction of the system with the solvent only. The theory is applicable to both low 
and high energy (with respect to solvent induced line broadening) vibrations. The results of this 
paper enable to qualitatively explain experimental observations on low energy vibrations presented 
in the preceding paper [A. Nemeth et al, arXiv:1003.4174vl] and to predict the time evolution 
of two-dimensional spectra in ultrafast ultrabroad band experiments on systems with high energy 
vibrations. 



I. INTRODUCTION 



Multi-dimensional ultrafast spectroscopic techniques 
open an increasingly broad experimental window into the 
dynamics of molecular systems on femtosecond and pi- 
cosecond time scales. Important experimental insights 
have been obtained by studying dynamics of small sys- 
tems such as vibrations of small molecules, peptides and 
proteins [l|, hydrogen bonding networks in water Q and 

geptides [3|, and electronic transitions in dye molecules 
1, |5j. At the same time, multi-dimensional spectroscopy 
have proven its strength on complex systems ranging 
from conjugated polymers ^6^, 't'], nanotubes f§\, protein 
chromophore complexes ISj-JJJ , all the way to solid state 
systems . The experimental work has been closely in- 
terlinked with theoretical efforts [l3-[13 which provided 
interpretation to this rich source of experimental infor- 
mation. One of the most interesting recent topics in 
this field was stimulated by the observation of oscilla- 
tory modulations of two-dimensional (2D) spectra, that 
were assigned to the presence and time-evolution of elec- 
tronic wavepackets [ill . [TtI [28j . While the oscillatory 
effects were predicted theoretically from exciton theory 
[13, Hoi, the details of the observed processes suggest that 
theories traditionally applied to model energy relaxation 
and transfer in these systems underestimate the role of 
coherence effects. For example, the decoupling of exci- 
ton population and coherence dynamics (secular approx- 
imation) seems not to be a valid approximation for the 
photosynthetic protein FMO [Tl| . 

Favorable properties of the studied multichromophoric 
systems, together with a partial success of exciton theory 
allowed to exclude vibrational coherence as a source of 



the oscillatory effects in the experiments on FMO 
In general however, a complementary understanding of 
the time-dependent spectral features of vibrational origin 
are clearly required in order to be able to disentangle 
the two effects. To this end, the topic of this and the 
preceding experimental paper, Ref. [2^ (paper I), is a 
system which can exhibit only vibrational coherence, a 
two-level electronic system interacting with well-defined 
intramolecular vibrational modes. 

For the description of non-linear spectroscopic signals 
of simple few electronic level systems with negligible re- 
laxation between the levels, the non-linear response func- 
tions derived in the second order cumulant approxima- 
tion (see e.g. (30| ) provide a framework of unprecedented 
utility. The non-linear optical response of such a system 
can be expressed in terms of line-broadening functions, 
which are related to the electronic energy gap correla- 
tion function [soj. Few electronic level systems provide 
models which enable to disentangle the content of exper- 
imental data otherwise difficult to interpret. For exam- 
ple, three pulse photon echo peakshift was identified by 
Cho et al. ^31] to reveal the energy gap correlation func- 
tion. Two-color photon echo peakshift has been found 
to measure excitonic coupling in a dimer fs?] or differ- 
ences between energy gap correlation functions of two 
chromophores [13 ■ Basic 2D spectral line-shapes were 
classified by Tokmakoff [s^l based on a two-level system, 
and finally, the relation between the non-linear spectro- 
scopic signals and the normalized energy gap correlation 
function was studied by de Boeij et al. [35] and for 2D 
spectra by Lazonder et al. (36j , to name but a few exam- 
ples. 

In some of these works, insightful information was ob- 
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tained by expanding the line-shape functions (or energy 
gap correlation functions) in terms of a suitable param- 
eter, e.g. the time delay between two interactions. In 
the present study, we follow the same general approach 
and expand the line-shape correlation functions (describ- 
ing the interaction of an electronic transition with an 
intra-molecular vibrational mode) in terms of another 
suitable parameter, the Huang-Rhys factor A/w. Here, 
A is the reorganization energy and w the frequency of the 
mode. The applicability of a low order expansion of the 
response functions in terms of the Huang-Rhys factor to 
the description of the 2D experiment is obviously system- 
dependent but for the dye molecule studied in Refs. jl^ 
and Q the qualitative results of the theory presented 
in this paper were confirmed by a full response function 
simulation of the 2D spectrum 

One of the first theoretical analyses of the vibrational 
line-shapes in electronic 2D spectra was done by Egorova 
et al. [M, 113 

using the density matrix propagation 
method. We complement this work by performing a de- 
tailed analysis of a two-level electronic system within the 
response function theory, where vibrational levels are not 
considered explicitely. Unlike in Refs. [2^ and [l^l we 
consider both rephasing and non-rephasing parts of the 
2D spectrum, and we study time dependent modulations 
of the 2D spectrum by slow vibrational modes. 

This paper is organized as follows. In the next section, 
we shortly review two-dimensional Fourier transformed 
spectroscopy and its theoretical description. In Section 
mil we introduce the non-linear response functions and 
the line-shape functions which form the basis of our de- 
scription of the 2D spectra. In Section ITVl we decompose 
the total 2D line-shape into several components accord- 
ing to their time dependence and derive elementary 2D 
line-shapes, which can be used to construct the total 2D 
spectrum. The time evolution of 2D spectra modulated 
by low and high energy modes are analyzed in Section [V] 
Our conclusions are summarized in Section [VII 



II. TWO-DIMENSIONAL FOURIER 
TRANSFORMED SPECTROSCOPY 

The principles of Two-Dimensional Fourier Trans- 
formed Spectroscopy (2D-FT) where described in detail 
in many references (see e.g. 0, [l5|) including the preced- 
ing paper, Ref. |29|. Here, we sum up results important 
for the theoretical modeling of the 2D-FT spectra in this 
paper. 

In 2D-FT we use three laser pulses with wave- vectors 
ki, k2, and ka arriving at the sample at times ri, T2, and 
T3, respectively. In case of sufficiently weak excitation 
pulses, the intensity of the signal observed in the direc- 
tion of — ki-|-k2-l-k3 is of third order in the total intensity 
of the incident light and the signal can be well described 
by third order time-dependent perturbation theory. The 
propagation of light through a material sample is influ- 
enced by the light-induced polarization, which is propor- 



tional to a material response function. The total third 
order response of the system to an excitation by three 
successive interactions with an external electric field of a 
laser is 

/ . X 3 4 

^ (^) E{^»(^3'*2,ii)-i?:(i3,i2,<i)}, (1) 

^ ^ n— 1 

where Q{t) is the Heaviside step- function. The Rn func- 
tions are defined according to the standard notation of 
Ref. [30], and represent different perturbative contribu- 
tions (know as Liouville pathways) to the total response. 
The response is a function of three time variables, where 
ti denotes the delay between the first and second interac- 
tion, t2 denotes the delay between the second and third 
interaction, and ^3 is the delay between the last inter- 
action and the measurement. The polarization of the 
sample at time ^3 after the arrival of the third pulse is 
given by a three-fold convolution of the response func- 
tion, Eq. JT]), with the incoming electric fields. For very 
short pulses (impulsive limit) one can identify times ti 
and t2 with the delays between the pulses (ti — |t2 — ti|, 
^2 = 1^3 — T2I) and one can approximate the polarization 
by those parts of the above expression that survive the 
so-called rotating wave approximation (RWA) (see e.g. 

Ref. 3). This yields Pi'^(i3,i2, ii) ~ 7^^'l(^3,^2,^l). 

(3) 

With different ordering of the pulses ki and k2, TZrwa 
contains different Liouville pathways. In particular 

P^'^Kh.h.h) w ~iR2{hMM) - iR^ihMM), (2) 

for Ti < T2 and 

n^^' (^3,^2, h) « -iRi{h,t2, h) - lRi{h, t2, h), (3) 

for Ti > T2. As pointed out in the Introduction, we are 
interested in two level electronic systems, for which the 
number of unique pathways is limited to four [30]. The 
electric field of the measured signal is proportional to the 
polarization as 

EfKhMM) ^ ^M^'KhMM). (4) 

n[ijj) 

In 2D-FT experiments, the signal travels into a direc- 
tion different from all incoming pulses, where it can be 
detected by a heterodyne detection scheme without in- 
terfering radiation from the three excitation pulses. The 
radiative prefactor Ld/n{u!) is removed and the 2D spec- 
trum is obtained by a double Fourier transform (experi- 
mentally only one of the Fourier transforms is performed 
numerically, the other is implicitly involved in the fre- 
quency resolved detection scheme) (1, Signals from 
both possible time orderings of the ki and k2 pulses (i.e. 
Ti < T2 and Ti > T2) are summed to form the 2D-FT 
spectrum yielding 

S{uJ3,t2..UJl) 
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+ FT^p [i?i(i3, t2,h) + Riih, t2, ii)](t^3, ^i), (5) 
where 

oo oo 

- J dt3 J dhf{h,h)e'^'''-'"''- (6) 



FTW[/(t3,tl)](^3,0.l) 

oo OO 

= Jdts I dhf{t3,h)e'"'''+'"''\ (7) 



The contributions stemming from different orderings of 
the pulses posses different phase factors with respect to 
times ti and ^3. Pathways R2 and R3 have a phase factor 
which turns equal to zero for all frequencies at ti = t^. 
They are referred to as rephasing pathways. Pathways 
i?i and R4 do not rephase in this sense and are referred to 
as non-rephasing. Thus, in the impulsive limit the 2D-FT 
spectrum can be calculated directly by applying a dou- 
ble Fourier transform to the rephasing and non-rephasing 
components of the third order response function. 

III. THIRD ORDER NON-LINEAR RESPONSE 
FUNCTION 

For a two-level electronic system, the third order re- 
sponse functions can be derived in terms of a single quan- 
tity called line-shape function, g{t). The line-shape func- 
tion is defined as 

t T 

g{t) =^ JdrJ dr'Cir'), (8) 


where C{t) is the so-called energy gap (or frequency- 
frequency) correlation function, describing the fluctua- 
tions of the electronic transition due to its interaction 
with the degrees of freedom (DOF) of the surrounding. 
Often it is advantageous to use another function, usually 
denoted as M{t), for the discussion of the influence of 
the bath DOF onto the spectrum (see e.g. [35|f36j). The 
M{t) function is (for high temperatures) equal to the real 
part of the correlation function C(t) normalized to 1 at 
t — 0. An analysis of nonlinear spectra based on M{t) 
was done also in the preceding paper, Ref. [2^. The 
response functions read in detail 

-R«(i3,t2,ti) = M|^e^"=«(*i-*^)+^"(*^'*^'*i), n = 2,3 (9) 



RnihMM) = |d|4e-*'"=«(*i+*-')+-^"(*^'*^'*i), n = 1,4, 

(10) 

where d is the transition dipole moment between the elec- 
tronic ground and excited state, and we defined 

fAhMM) = -g{ti)~g*{h)~9*{t2)+g*{t2+h) 
+ gih + t2)-giti+t2 + h), (11) 

f2{t3, t2, h) = ~g*{h) - g*{h) + g{t2) - g{t2 + h) 

~g*iti + t2)+g*ih+t2 + h), (12) 

hihMM) = -g*{ti)-g{h)+9*{t2)-g*{t2+h) 
-g*{h + t2)+g*{h+t2 + h), (13) 

/4(i3, t2,h) = -g{h) - gits) - g{t2) + g{t2 + is) 

+ .9(^1+^2) -3(^1+^2+^3). (14) 

Using Eqs. ^ to we can calculate the impulsive 
response of a two-level electronic system interacting with 
a rather general environment that is described by some 
correlation function C{t). 

One of the most common models for the correla- 
tion function C{t) is the so called Brownian oscillator 
model, which enables to describe both overdamped sol- 
vent modes as well as underdamped intramolecular vi- 
brational modes. The Brownian oscillator model and its 
application to calculation of optical spectra is extensively 
described in Ref. (soj . In our considerations, we repre- 
sent certain interesting intramolecular modes by the un- 
derdamped version of this model, assuming further that 
the damping of the oscillator is negligible. For the de- 
scription of the experimental results, Ref. [2^, it seems 
to be a good approximation for times less than 1 ps. The 
corresponding energy gap correlation function thus reads 

M 

C{t) — AajS(T) coscjt — lAwsinwi, (15) 

where S(T) = coth^^^l^)- ^'^^^^ ^q. ^ it follows 
that 

g{t) = iXt + i- smujt + -E{T) [1 - coscji] . (16) 

Here, A is the reorganization energy and u is the fre- 
quency of the mode, respectively. The overdamped Brow- 
nian oscillator model is used for modeling the solvent cor- 
relation function. We will be interested in the properties 
of the corresponding solvent (7(i)-function for t > Tc only, 
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where Tc is some characteristic correlation time of the sol- 
vent. Particularly, if the solvent energy gap correlation 
function does not have any slow static component, i.e. 
if so-called static disorder of the system is negligible, we 
can write for t > Tc (see Ref. ^30j ) 



-I- Gn{ti,t3)s'mujt2, 



(22) 



git) = gconst + g{t), 

where g{t) has the property 

git + T) =git) + gir), 



(17) 



(18) 



i.e. for sufficiently long times, the function g{t) is linear. 
This property will be use below to simplify Eqs. (fTTj) to 



IV. MODULATION OF NON-LINEAR 
RESPONSE BY VIBRATIONAL MODES 



Third order response of weakly coupled 
vibrational modes 



Let us consider an electronic two-level system that is 
coupled to an environment consisting of one dominant 
underdamped vibrational mode and a bath with a macro- 
scopic number of DOF. The latter DOF are represented 
by an overdamped Brownian oscillator model and its 
corresponding line-shape function gbath{t) which has the 
property of Eq. ([T7]). The total energy gap correlation 
function of the system will be assumed in the form 



g{t) = gvtb{t) +gbath{t), 



(19) 



where gvib{t) is defined by Eq. ([TC)) . The response func- 
tions, Eqs. ^ and ^TU\\ . can thus be factorized into pure 
solvent and pure vibrational parts. Since we deal here 
with one electronic transition only, we assume d — I for 
the rest of the paper. If the Huang- Rhys factor X/u) (i.e. 
the coupling of the vibrational mode to the electronic 
transition) is sufficiently small, we can expand the vibra- 
tional contribution and write to first order in X/uj 



1 + -F„{h,t2,t3] 



2,3 



(20) 



with K , H, and G functions discussed below. We used 
the fact that if Eq. ([TBI) inserted into Eqs. ([11]) to 
(fTil) . the first two terms we obtain are fn{t3,t2,ti) — 
iX{t3 ± ti) — 2AS(r)/tj -I- . . .. These two terms are kept 
in the exponential of Eqs. ([20]) and pi|) . and the rest is 
expanded in terms of X/uj. The only formal differences 
between the rephasing and non-rephasing response func- 
tions, Eqs. ([20]) and ([2T|) . are the different phase factors 
in times ti and ^3. 

The expansion of the exponential function in Eqs. ([20)) 



and ([21]), exp{^F„} I + Ai^„, requires ^F^ to be a 
small number. The function Fn is an oscillating function 
bound between certain values which represent the points 
of the extreme deviation from the exact expression for 
the response functions. One of such values occurs e.g. at 
ti = and t3 = 0, where F„ = 2E{T). The validity of the 
approximation is therefore limited also by temperature. 
For high temperatures, however, i.e. \i hw ^ 2kBT and 
correspondingly S(T) « 1, the validity is limited only 
by the value of the Huang-Rhys factors. The 2D spectra 
that are the main subject of this paper are obtained by 
Fourier transforms in times ti and t^. In applying the 
approximation it is therefore important that the approx- 
imated response function retains the correct periodicity 
in ti and is . The expansion can be generalized for several 
vibrational modes, but one has to bear in mind that the 
factor limiting the validity of the approximation might be 
the sum of Huang-Rhys factors of all expanded modes. 
Using the well-known goniometric relations 

sin(a -I- /3) = cosasin/3 -|- sinacos/3, (23) 



cos(a + /?) = cosacos/? — sinasin/3 (24) 

the if, H , and G functions can be expanded into a lin- 
ear combination of simple cos and sin functions. The 
actual forms of these functions can be found in Ref. [H 
or in Appendix [A) The t2— dependence of the response 
related to vibrational modes can thus be approximately 
separated into stationary, sine oscillating and cosine oscil- 
lating parts, using relatively simple factors of Eqs. ([Aip 
to (|A12[) from Appendix [Aj 

Due to the properties of the solvent line-shape function 
gbath{t), namely Eqs. ([T71) and ([T5)) . the response of the 
solvent is i2-iiidependent for long times t2. Thus, if t2 > 
Tc we can write 



l + -Fn{tut2M) 



n = 1,4. 



(21) 



Here, the i?n(t3, ^2, ^i) functions comprise the pure sol- 
vent response (except of the optical frequency factor 
g-i<^e,(t3±ti)^ cf. Eq. ([IDanddn])), and 

F„{ti,t2,h) = Kn{tl,h) + Hn{tl,h) C0SUjt2 



-gbath{ti)-gbath{t3) 



(25) 



MhMM) = MhMM) = e-5-".(*i)-9--(*3). (26) 

This enables us to drop the i2-dependence of the solvent 
response function in the following discussion, bearing in 
mind that we always assume t2 > Tc- 
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B. Two-dimensional Fourier-transformed spectra 

The 2D-FT spectrum was defined in Section |lll By tlie 
double Fourier-transforms, Eqs. ([5]) and ([7]), we obtain 
the total spectrum and its rephasing and non-rephasing 
parts at a given delay t2- The expansions in Eqs. (PO)) 
and (PT|) suggest that the 2D spectrum is composed of a 
vibration-free spectral shape originating from the solvent 
response function R and some vibrational contributions. 
The pure solvent contributions read 

xe^"-('^-*^)](c.3,Wi), n = 2,3, (27) 



Now we will use elementary rephasing and non- 
rephasing 2D spectral line-shapes, Eqs. and ([5T|) . 
to construct a dependent line-shape of a molecule 
possessing a vibrational mode with a sufficiently small 
Huang-Rhys factor. To this end, we use the fact that 
the vibrational part of the response function consists of 
the t2— independent functions Kn{ti,t3), Hniti^t^) and 
Gnitijt^) that contain only sines and cosines of vari- 
ous combinations of arguments ojti and tuts. The whole 
^2— dependence is contained in Eq. ([2^ . While perform- 
ing a Fourier transform of the response functions, Eqs. 

and (I^T|) , we Fourier transform the solvent response 
function multiplied by the cosine and sine terms from 
the functions Kn{ti,t^), i?„(ii,<3) and Gn[ti,ti)- The 
total rephasing spectrum can thus be constructed from 
six possible vibrational contributions. We have e. g. 

Z(3)(u;3,C^l) = J7T2(^)[e»("=.+A)(ti-t3) 



X e-*"='(*i+*^)](a;3,cji), n = l,4. (28) 
The solvent 2D-FT spectrum is therefore given by 
S'soii>(w3,t^i) = 2'2°'"(a;3,a;i) -|- Z3°'"(a;3, wi) 

-f yr'"(^3,^i) + i^r"'(c^3,^i). (29) 

Considering Eqs. and (pS)) we see that one can write 



Z(a;3,a;i) 



Zr'"(c.3,wi) = a(c^3)e*(wi), 



y(c^3,^i) = ll'''"'(c^3,wi) 



rr"'(^3,^i) = 5(^3)5(^1), 



where we defined 



oo 

Q{u]) = j dt e' 



(30) 



(31) 



(32) 



Solvent contributions to the 2D spectrum can thus 
be expressed using just two elementary 2D line-shapes 
Z[uj3,uji) and y(a;3,a;i), which can be further simpli- 
fied. The real part of Eq. is directly proportional 
to the absorption at frequency w. Consequently, one 
can conclude that the real part of the total 2D line- 
shape, Re(Z + y) = Re(C/(w3))Re(^/(cL>i)), represents an 
absorption-absorption correlation plot, while its imag- 
inary part, Im(Z -\- Y) — Im(^/(w3))Re(^?(a;i)), corre- 
sponds to an absorption-refraction correlation plot. For 
simple symmetric line shapes with a maximum at some 
frequency iOmaxi e.g. for a Lorentzian, Z and Y are re- 
lated by flipping of the line shape with respect to the 

Wl-axis at UJi = UJmax- 



X i?2 (^3 , i 1 ) COS Ujt^i] {UJZ , Wi ) . 



(33) 



The remaining expressions of the same kind are summa- 
rized in the first part of Appendix |B] From now on, 
we use only ^2(^3,^1), because R^itziti) = Rsihjti) for 
^2 > Tc (cf. Eq. (|26p ). and we drop the lower index 
for the line-shapes Z and Y. Similarly we obtain shapes 
Y^os, Yco} etc. for the non-rephasing part by applying 
FT^P to the Ri R4) function. 

The vibrational frequency uj is much smaller than the 
frequency of the visible excitation light at which the 2D 
spectrum is measured and consequently the sine and co- 
sine in the Fourier transform introduce only shifts of the 
spectral feature without significant changes of its shape. 
Using the relations 



cos ujt 



and 



we obtain e.g. 



^cos(w3,wi) = -^^(^3 - A - - A) 



sinwi — i- 



(34) 



(35) 



+ 2^('^3 - A + - A), 



(36) 



— Z(a;3 — A + cj, wi 



A), 



(37) 



and analogously for other shapes. All shapes can thus be 
written as simple linear combinations of the shifted basic 
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shapes Z and Y . To save writing out of all arguments, 
we define 



and 



fis = - A, ill ^ - A, 



(38) 



and the cosine oscillating rephasing part composed of the 
Fourier transforms of the i?2(^3,^i) and H^{t^,ti) func- 
tions can be written as 



A 



Z(3) + ^(1) _ ^(3+1) _ z 



UJ 



'{-+) = 



(39) 



(40) 



etc., where upper indices +, — and of the newly defined 
Z functions indicate the additional shift by frequency 
-foj, —UJ or X oj, respectively. Using the above defini- 
tions we can write all the contributions of the actual 2D 
spectrum e. g. as 

zill{n„n,) = l{z(-^^ + z(+^^ 



+ i- 



.(1) 

"'sin 



.(3+1) 

sin 



(45) 



The non-rephasing contributions are similarly based on 
the non-rephasing functions with subindices 1 and 4 and 
read 



Y - i-Y 



(1) 



, . '(-^) (.^cos ^~ ^cos ) 



(46) 



zl^)(173,f^i)-^(^(-"'-^(+"^). (41) 

Again, the remaining equations are summarized in Ap- 
pendix Ib] 

According to Eqs. (^0)1 to ([^ (and assuming t2 > t^) 
the total response can be split into a independent 
part (solvent and vibrational contribution given by func- 
tions Kn{h^ti) defined in Eqs. ((1T|) to ((XH ) and a 
part oscillating with t2 (given by functions Hn{tz,ti) 
and Gn{h,ti) of Eqs. to (|A12p ). Correspondingly, 

the total 2D-FT spectrum can then be split into a t2- 
independent part, 5*', and parts, 5^°*, and S**™, oscillat- 
ing with cosa;i2 and sincji2, respectively, 

5(113, t2,^l) = 5°(ll3, f^l) + 5™^(ri3, f^l) C0SUjt2 



-hS'""(fl3,f^i)sina;t2- (42) 

Collecting all parts of the response, we can express all 
contributions to the total 2D-FT spectrum in terms of 
the Z and Y line-shapes defined in Eqs. ((30)) and (|3T|| . 
We keep the lower indices R and NR to denote the 
rephasing and non-rephasing parts of the spectrum, re- 
spectively. The static rephasing part consisting of the 
Fourier transform of the rephasing solvent response and 
functions ^^2(^3,^1) and -^^3(^3,^1) reads 

Sl{n,^n,)^\z + ^-Z^2^ + -^{T){Z[]^, + Zf^± (43) 

The sine oscillating rephasing part which is based on the 
^2(^3,^1) and 63(^3,^1) functions reads 



5r(f^3,f^i)«-- 

LiJ 



yd) , ^(1) _ ^(3+1) 
sin ' sin sin 



.A r 



ZW _ ;7(3+i) 



(44) 



5^^(f73,f^i)«-S(r) 

UJ 



y(l) , y(3) _ y(3+l) 

sin ' sin sm 



A 

■ I— 



Y(1) _ y(3-Hl) 



(47) 



5^1(173,^^1) «—s(r) 



y(3) _Ly(i) _y(3+i) _y 



.A r 



yd) _ Y 

sin sin 



(3+1) 



(48) 



The above equations enable to characterize completely 
the total time-dependent 2D-FT spectrum of a two-level 
electronic system weakly coupled to a vibrational mode. 



V. ANALYSIS OF TIME-DEPENDENT 2D 
LINE-SHAPES 

In this section, all calculations are performed assum- 
ing T = 300 K and the bath correlation function in the 
form of an overdamped Brownian oscillator with reorga- 
nization energy Xhath = 100 cm~^ and correlation time 
Tbath = 150 fs. These numbers represent some character- 
istic values motivated by the molecular system studied in 
Paper I and Ref. fB^. For the Brownian oscillator model, 
the long time line-shape function has a constant term 
dconst — —XTc{2kBTTc/h — i) and a linear term, which 
can be written as g{t) = —gconstt/Tc- The high temper- 
ature limit ksT ^ h/Tc is valid for these parameters. 

If the influence of vibrational modes converges to zero, 
i.e. for A/w — >■ 0, the total spectrum is formed by the 
sum of the non-shifted rephasing shape Z and the non- 
shifted non-rephasing shape Y . Qualitatively, these line- 
shapes are the same as those presented in Fig. [T] The 
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Figure 1: The stationary contributions (from top down) 
S%{u!3,uji), S%ii{ijJ3,uji) and their sum S^{oj[i,u!i) for the 
modulation by a low energy vibrational mode of cj = 140 
cm~^ and A = 80 cm~^. The figures on the left present the 
real and the figures on the right the imaginary parts of the 
spectra, respectively. Apart from a slight asymmetry of the 
line-shapes with respect to the diagonal or anti-diagonal line, 
these figures represent very well the pure solvent line-shapes, 
too. 



real rephasing part of the spectrum has a characteris- 
tic orientation of the positive peak along the diagonal 
line, while the real non-rephasing part is elongated along 
the anti-diagonal line. Both line-shapes have negative 
"wings" stretching perpendicularly to the axis of elonga- 
tion of the particular spectrum. If combined, they form 
a symmetric shape of the total response with solvent in- 
duced broadening only. The dispersive (imaginary) part 
of the solvent spectrum is also formed by the rephasing 
and non-rephasing elementary 2D line-shapes which are 
characterized by a nodal line parallel to the elongation 
axis of their real counterparts. Thus, the rephasing and 
non-rephasing dispersive parts differ from each other by 
the orientation of their nodal lines. The combined spec- 
trum has a nodal line parallel to the wi-axis. 

The 2D line-shapes Z and Y represent the basic forms 
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Figure 2: Demonstration of the construction of the 2D line- 
shape S'lj"' (oJs , wi ) . Real parts of the individual contributions 
are depicted in panels (a) to (e), together with the signs of the 
particular contributions. The total line-shape 5^°" (ctJs , ctJi ) is 
depicted in panel (f). The central negative feature originates 
from the unshifted line shape Z(Loz-,io\). The signs in the 
lower right corner in panels (a) to (e) represent the signs with 
which the diagrams contribute in Eq. (|45|) . The numbers in 
panel (f) denote the predicted relative intensity of the fea- 
tures. 



from which all spectral components can be constructed 
by means of a shift by the vibrational frequency cj. We 
demonstrate this on the contribution 5'^°* (cf. Eq. (|45p ). 
Fig. [2^-e depicts the individual terms in Eq. P5)) ac- 
cording to Eqs. ([41]), ((B7| and (|B9|) . and^f) presents 
the total contribution 5"^°^. The oval shapes represent a 
single contour of the real part of the rephasing line-shape 
Z, thus showing their characteristic orientation and posi- 
tion in the spectrum. The numbers Fig. [2f represent the 
intensity of the corresponding spectral feature. Empty 
contours represent negative peaks, whereas shaded con- 
tours indicate positive peaks. From this schematic pic- 
ture one can draw several conclusions. First, if the shift 
uj is small and all contributions overlap significantly, the 
positive part of the sum is elongated along the diagonal, 
with asymmetric negative parts along the anti-diagonal 
line. The negative features are stronger in the lower right 
corner of the 2D spectrum. Second, for vanishing w the 
total contribution tends to zero, because all terms cancel. 
Consequently, for the same Huang-Rhys factor, the mod- 
ulation by slower (lower energy) modes is weaker than 
the corresponding modulation by faster (higher energy) 
modes. Third, for uj larger than the width of the ele- 
mentary line-shapes Z, the rephasing cos contribution 
has the strongest positive feature on the right from the 
central frequency, a strong negative feature in the lower 
right corner, while in the center the contribution is neg- 
ative. In the following sections we will demonstrate this 
prediction by numerical evaluation of Eqs. to 
Similar intuitive pictures can be created also for the non- 
rephasing spectra, using the pictorial representation of 
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the line-shape Y. 

We will distinguish between fast (high energy) and 
slow (low energy) vibrational modes relative to the in- 
verse homogeneous spectral width Agoiv (FWHM) of the 
solvent contribution. If a; Asoiv, i-e. vibrational cross- 
peaks in the 2D-FT spectrum or vibrational lines in the 
absorption spectrum would be indistinguishable, we as- 
sume the vibrational mode to be slow. The opposite case, 
LJ S> Aso/.y. characterizes fast vibrational modes for which 
individual vibrational features could be distinguished in 
optical spectra. 




(D ffi (D 



CO] /m ODi /m 

Figure 3: The stationary contributions (from top down) 
^^(cja, cji), S%ji{uJ3,uJi) and their sum S'^{l03,lji) for the 
modulation by a high energy vibrational mode of cj = 1000 
cm" ^ and A = 300 cm~^. The figures on the left present the 
real and the figures on the right the imaginary parts of the 
spectra, respectively. 

Although the vibrational contribution is mostly dy- 
namic, i.e. dependent, there is a certain time- 
independent contribution to the total 2D spectrum that 
originates from the vibrations. The stationary parts of 
the 2D spectra for modulation by a low energy vibra- 
tional mode of a; = 140 cm"^ and a high energy vibra- 
tional mode of a; = 1000 fs are presented in Figs. [T] and 



131 respectively. The reorganization energies have been 
set to A = 80 cm~^ and A — 300 cm"'^, respectively. In 
case of the low energy mode, the shift of the line-shape 
by w does not lead to any crosspeaks separated from the 
main diagonal peaks. The observed line-shape is only 
slightly broadened and distorted from the pure solvent 
line-shape by the vibrational contribution. In case of 
the high energy mode, we can distinguish separate vi- 
brational cross- and diagonal peaks. In all cases we can 
identify the characteristic direction of elongation of the 
2D line-shapes. The imaginary parts of all spectra also 
show characteristic orientations of the nodal lines. We 
can conclude that, in the static part of the 2D spectrum 
modulated by a high energy mode, the line-shapes are 
dictated by the chromophore solvent interaction, while 
the positions of the peaks are given by the vibrational 
frequency. 



A. Modulation of 2D spectra by a low energy mode 

The time evolution of the 2D spectra is mediated 
by the sin(u;<2) and cos(wt2) functions which vary the 
magnitude of the contributions in Eqs. (gl]), (gSJ, (gT]) 
and dMl). A closer look at Eqs. to (gS]) reveals 

that the rephasing and non-rephasing contributions have 
a very similar functional form. They differ mainly by the 
exchange Z < — > Y. Apart from this exchange, compar- 
ing the two sin contributions, Eqs. (|44| and (|47)) . and 
the two cos contributions, Eqs. (|45p and pS)) . there is 
a flip of the sign in their first terms. This suggests that 
the rephasing and non-rephasing modulations of the to- 
tal 2D shape differ not only in the orientation of their 
elongation, but also in phase. For cosa;i2 = 1, the to- 
tal rephasing contribution is at its maximum, while at 
the same time, the non-rephasing one is at its minimum. 
Due to the third term in Eqs. (gl]), (gS]), (gT]) and (gHD 
this property is not exact. However, by studying dia- 
grams similar to Fig. [2] for the various contributions, we 
have concluded that a cos non-rephasing line-shape can 
be obtained from the rephasing one by turning the pic- 
ture by tt/2 anti-clockwise (corresponding to Z i — > Y 
exchange), mirroring it by the diagonal line, and chang- 
ing the sign. 

Figs, g] and [5] present the cos and sin oscillating parts 
of the total 2D spectrum for modulation by a low energy 
vibrational mode. We can immediately notice that for 
both the sin and the cos modulation the real rephasing 
and non-rephasing contributions are approximately out 
of phase. For the set of parameters used here, the maxi- 
mum amplitude of the cos contribution is approximately 
twice the maximum amplitude of the sin contribution, 
so that the cos oscillation is the most significant fea- 
ture in the i2-evolution of the 2D spectrum. Neglecting 
therefore the sin contribution for a while, we can pre- 
dict a mutual out of phase modulation of the rephasing 
and non-rephasing parts of the 2D spectrum. This leads 
to the characteristic oscillation of the diagonal and anti- 
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Figure 4: The cos contribution for the modulation by a low 
energy vibrational mode of oj = 140 cm~^. From top down: 
the rephasing and non-rephasing parts, and the sum spec- 
trum. The figures on the left present the real and the figures 
on the right the imaginary parts of the spectra, respectively. 
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Figure 5: The sin contribution for the modulation by a low 
energy vibrational mode of a; = 140 cm~^. From top down: 
the rephasing and non-rephasing parts, and the sum spec- 
trum. The figures on the left present the real and the figures 
on the right the imaginary parts of the spectra, respectively. 



diagonal width of the total 2D spectrum (of. Fig. |6]), 
which was demonstrated experimentally in the preceding 
paper, Ref. [1^. The nodal line rotation in the disper- 
sive part of the total 2D spectrum has a similar origin. 
As the rephasing and non-rephasing parts of the disper- 
sive contribution have different orientation of their nodal 
lines, the total 2D spectrum exhibits oscillations of the 
nodal line angle, due to the change of relative amplitudes 
of rephasing and non-rephasing parts. 

Fig. [S] demonstrates how the combined non-rephasing 
and rephasing signal gives rise to the characteristic ex- 
perimentally observed modulation of the 2D line-shape. 
The results of our approximate analysis in this section are 
well confirmed by the experimental 2D spectra presented 
in Ref. (29| . Quantitative agreement with the experi- 
ment was achieved in Ref. [5] by modeling with the full 
expression for the corresponding response function, Eqs. 
® to HID. 



B. Modulation of 2D spectra by a high energy 
mode 

Encouraged by the successful comparison of the the- 
ory presented here with an experiment involving low en- 
ergy vibrational modes, we can extend the treatment to 
high energy modes. For fast vibrational modes, where 
huj ^ ^ksT, the validity of the expansion in Eqs. (PO)) 
and is limited only by the value of the Huang-Rhys 
factor A/cj, because coth(fta;/2fcBT) « 1. The composi- 
tion of different contributions to the 2D spectrum in this 
case, can be easily understood directly from considera- 
tions similar to those which led to Fig. [21 The rephasing 
real part of the cos contribution shown in Fig. [7] corre- 
sponds very well to the depiction of Fig. [5f. A negative 
peak is present in the center of the spectrum, while two 
positive peaks are below the diagonal and one above the 
diagonal. The stronger of the two positive peaks below 
the diagonal is positioned at higher uji frequencies. All 
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Figure 6: The time evolution of the total 2D spectrum of an 
electronic transition modulated by a low energy vibrational 
mode of ct; = 140 cm~^. From top down: 2D spectra at times 
T = 0,^, ^, and fs. The fi gures on the left present the 
real and the figures on the right the imaginary parts of the 
spectra, respectively. 



other contributions can be constructed in a similar way, 
including those oscillating with sin{ujt2) (cf. Fig. [5]). 

Figure [H] presents the time evolution of the total 2D 
spectrum at different phases of the oscillation. An in- 
teresting observation is that for high energy vibrational 
modes the diagonal peaks oscillate with the opposite 
phases than in case of the low energy mode modulation, 
and with the opposite phase with respect to some of the 
cross-peaks. The shape of the cos contributions for low 




(D 

coi /m * 



Figure 7: The cos contribution for the modulation by a high 
energy vibrational mode of a; = 1000 cm^^. From top down: 
the rephasing and non-rephasing parts, and the sum spec- 
trum. The figures on the left present the real and the figures 
on the right the imaginary parts of the spectra, respectively. 



and high energy mode modulations (Fig. |4]vs. Fig. [7]), 
which are constructed from the same Z and Y but with 
different amount of shift from the transition frequency, 
again enables to understand the difference. Unlike the 
central diagonal peak, the crosspeaks oscillate between 
positive and negative values. The phase of these oscil- 
lations enables, for example, to clearly distinguish the 
2D spectra in times <2 = 7r/2aj and t2 — 37r/2aj by the 
real part of the spectrum only, without referring to the 
imaginary parts. For low energy modes, these two phases 
might be very difficult to distinguish. 

The presence of negative features in different parts of 
the 2D-FT spectrum might influence the measured 2D- 
FT spectrum even if the spectrum of the pulse is not 
able to cover vibrational crosspeaks (see e.g. Ref. [37|). 
However, sufficient information for the characterization 
of wavepacket motion cannot be obtained from such a 
measurement. If the experiment is performed with a laser 
pulse of a bandwidth comparable with the width of the 
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Figure 8: The sin contribution for the modulation by a higli 
energy vibrational mode of oj = 1000 cm~^. From top down: 
the rephasing and non-rephasing parts, and the sum spec- 
trum. The figures on the left present the real and the figures 
on the right the imaginary parts of the spectra, respectively. 



vibrational peaks in the absorption or 2D spectrum (such 
as in Refs. [1, [l^ and (s^l) the frequency of the laser 
pulse relative to the main diagonal peak determines the 
amount of influence the fast mode has on the 2D spec- 
trum. In particular, in Refs. [5L[2^ the central laser pulse 
frequency is lower than the position of the main diagonal 
peak, and the 2D spectrum is thus measured away from 
the fast vibrational features. Any vibrational modes with 
vibrational frequency significantly faster than the width 
of the laser pulse can thus be ignored from the evalu- 
ation, which was successfully done in Ref. [5i]. On the 
other hand, if the laser pulse falls into the spectral region 
between the diagonal peaks, as is the case in Ref. |37| . 
the 2D spectrum is distorted by the high energy mode 
even if its frequency is larger that the width of the laser 
pulse. 

Application of ultra broadband pulses that would cover 
significantly the whole region where a studied vibrational 
mode absorbs, promises to enable a complete character- 



ization of the photo-induced vibrational wavepacket mo- 
tion. Our results present reference spectra for an un- 
derdamped wavepacket motion in a harmonic potential. 
Distortions from this reference case can reveal changes in 
the shape of the wavepacket due to weak relaxation, an- 
harmonicity, and also differences between the ground and 
excited state potential energy surfaces. From this point 
of view, 2D-FT experiments are of interest not only for 
molecules in solutions, but also for gas phase molecules, 
were they can characterize fundamental properties of vi- 
brational DOF. 



VI. CONCLUSIONS 

We have presented a qualitative theory of the time 
dependence of two-dimensional line-shapes of a two- 
level electronic system weakly coupled to a vibrational 
mode. We have studied the time evolution of the two- 
dimensional spectra in case of low and high energy vi- 
brational modes and showed that the two-dimensional 
spectrum at t2-times larger than the correlation time of 
the solvent can be qualitatively predicted based on the 
purely solvent two-dimensional line-shapes and knowl- 
edge of the vibrational frequency and Huang-Rhys factor. 
The theory predicts oscillations of the diagonal and anti- 
diagonal width of the two-dimensional spectrum caused 
by the rephasing and non-rephasing signals being modu- 
lated with opposite phase. In case of a low energy vibra- 
tional mode, the signal is predominantly modulated by 
a cosine term, as confirmed in the accompanying exper- 
imental paper. For high energy modes, where the line- 
shape width is smaller than the vibrational frequency, 
our theory predicts cross-peaks oscillating similarly to 
the low energy mode case, but the diagonal peaks oscil- 
lating with a phase opposite to the cross-peaks. Ultra 
broadband laser pulses would be necessary to observe 
high energy modes of molecules in liquid phases. How- 
ever, in gas phase, a full characterization of vibrational 
motion based on two-dimensional spectroscopy would be 
possible. 
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Figure 9: The time evolution of tlie 2D spectrum of an 
electronic transition modulated by a high energy vibrational 
mode ofuj = 1000 cm~^. From top down: 2D spectra at times 
T = 0, ^, ^, and fs. The fi gures on the left present the 
real and the figures on the right the imaginary parts of the 
spectra, respectively. 



Appendix A: Components of the vibrational 
response 

In this appendix, we present explicite formulae for the 
components of the vibrational response, Eq. (P^ . The 
t2-independent part Kn{ti,t^) of Eq. (|22D reads 

Ki{ti,tj,) — ~i[smLjti — sin wis] 



+ S(T)[coswti + coswts], (Al) 
^2(^1,^3) = i[sinujti + sincjia] 

+ S(T) [cos uti + cos ujtz] , (A2) 
-^^3(^11^3) ~ i[sinu)ti — sin wis] 

+ S(T) [cos uti + cos utz] , (A3) 
K4^{ti^t^) = —i[smLuti +sincjt3] 

+ E{T)[cosujti + cos wis]. (A4) 

The cosine oscillating part i/„(ti, ^3) cosa;t2 of Eq. 
can be written using 

Hi{ti,tj,) — S(r)[l — cos — coswii + cosa;(ii + t^)] 
+ i[sinwti — sincjis — sinw(<i + is)], (A5) 

H2{ti, ts) = — S(r)[l — cos w<3 — coswii + cos w(ti + t^)] 
+ i[sinwti — sincjis — sinw(<i + t^)], (A6) 

Hsitijt^) = — S(r)[l — COS w<3 — cosujti + cos w(ti + t^)] 
+ i[sinwti + sinwis — sinci;(fi + t^)], (A7) 

H4{ti,t3) — E{T)[1 — cosujts — cosujti + cosa;(ii + t^)] 

+ i[siiiujti + s'mujt^ ~ sina;(ii + ^3)]. (A8) 

Finally, the sine oscillating part G„(ti, is) sina;i2 of Eq. 
has a form of 

Gi(ii, ts) — S(r)[sinwt3 + sincjii — sinw(<i + t^)] 

+ i[l ~ coswts + coscjii — cosa;(ii + t^)], (A9) 
^2(^1,^3) = -S(T)[sinwi3 + sinwii - sinw(ti + t^)] 

+ i[l — cosLdt^ + coswti — coscLi(ti + is)], (AlO) 
G'3(ii, ts) — — S(T)[sinwt3 + sinwii — sinw(ti + t^)] 

— i[l — coswis — coswti + cosw(ti + is)], (All) 
G4(ii, is) = S(r)[sin wis + sinwii — sinci;(ii + is)] 

— i[l — coswis — coswii + cosa;(ii + is)]. (A12) 
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Appendix B: 2D line-shape components 

This appendix lists some intermediate definitions used 
to derive Eqs. (gS]) to (gS]) for the time dependent 2D 
hne-shapes. Line-shape components complementing Eq. 
read 



X R2{t3,ti) C0SUjti]{uJ3,UJi), 



(Bl) 



l-i(^^3,f^i) = ^(^("+^-^(°-)), (B7) 



Z(L+^)(f]3,ni) = i(z(+-'+^(-+)), (B8) 



ziir'\ns,n,)='-{zi-+^-Zi+-^), (B9) 



(BIO) 



(B2) 



(Bll) 



X i?2(i3,ii)sina;ii](cj3,iui), 



(B3) 



X i?2(i3,ii)cosa;(t3 +ii)](w3,cJi), (B4) 



(B12) 



(B13) 



X R2{t3,ti) sinw(t3 + ti)]{uj3,ui). 



(B5) 



Lineshape components complementing Eq. (|4ip read 
zW(l]3,f^i) = ^(^("+)+^("-)), (B6) 



y£+i)(f^3,f^i) = I , (bw) 



(3+ 



y( — ) _ y(++) 



(B15) 
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